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In this work we present a theoretical study of EELS (electron-energy-loss spectroscopy) ex- 
periments on the Ceo molecule. Our treatment of the problem is based on the simple two-fluid 
modeJ-|Originally proposed for the description of plasma oscillations in graphite and fullerenes (see 
RefsOtfl). It is shown that in spite of the simplicity of the model the calculated intensities of the 
EELS peaks are in good agreement with experimental data which may indicate that the model can 
be used as a simple and effective tool for the investigation of the collective behaviour of electrons in 
fullerene systems. 
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I. INTRODUCTION 

Investigations of plasmons in small clusters of nanometer size range have a long history. Originally plasma oscilla- 
tions were studied in small (R = 2 10 nm) metal ^articles. It turned out that plasmons in metal clusters can be well 
described by linearized hydrodynamic equationsETcl. The positions of plasmon peaks, their number and fast electron 
cross-sections calculated within these equations are in good agreement with experimental data. New possibilities of 
the investigation of plasmons in small objects appeared after the discmrery of fullerenes and the cheap methods of their 
synthesis n Plasmon oscillations were observed in fullerene moleculesucl, carbon nanometer size tubesEI and multishell 
fullerenesB. 

An important feature of fullerenes is the fact that these molecules are more or less of a spherical form and their 
electrons can be considered to be confined to the surface of the sphereO. Hence, it seems to be a good starting point 
to consider the Cgo molecule as a rolled up graphite plane. This approximation is expected to work well for the 
description of collective excitations. 

The graphite plane consists of regular hexagons formed by carbon atoms. Each carbon atom is connected by 3 a- 
bonds with its three neighbours, meanwhile 7r-electrons form a common 7r-system consisiting of orbitals perpendicular 
to the plane. In linear response theory it and tr-electron subsystems can be considered to be independent with respect 
to the electrical field parallel to the graphite plane because dipole transitions from 7r to cr-orbitals are prohibited. 
This means that we can describe the plasmons in graphite whose wave vector is parallel to the plane as the plasmons 
in a "two-liquid" electron system; it is proposed that both the liquids move in a common average potential. 

The graphite loss function -lm(l/&i,(|Ui)) with respect to the external electric field perpendicular to the c-axis shows 
two peaks at 7 and 28 eV (see RefsJlMTl3). They are called n- and a- (or "7r + er"-) plasmons. Strictly speaking such 
a division is rather conventional. In fact, both n- and tr-electrons make a nonzero contribution to er- and 7r-plasmons 
respectively. i— ir-i 

In the articlesOcJ Cazaux used the two-fluid model in order to explain the results of EELS studies on graphite. 
The same idea was used in Ref.E£l for the description of plasmons in the Ceo molecule: it was proposed that 7r-electrons 
be affected only by induced electric potential, meanwhile cr-electrons be governed also by a restoring force .-This force 
appears in the model via an effective description of a contribution of valent electrons to the polarizabilityEOiij. 

Using the two-fluid model Barton and Eberlein predicted the existence of plasma oscillations in the Ceo molecule 
which can be classified by the spherical quantum numbers I and m and the index j (j = 1, 2) which corresponds to 
7T- and CT-plasmonal3. However, the_alasmons with large I cannot exist because of their strong damping. 

EELS studies on Ceo (see Refs.cKl) show a peak at 6.5-7 eV (also called "7r-plasmon" ) and a rather broad "er- 
plasmon" peak at 15-16 eV . The former peak is identified usually as the I — 1 7r-plasmon, meanwhile the latter 
peak is considered to be a combination of the I = 1 er-plasmon and probably the plasmons with I > 1. However the 
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existence of multipole plasmons in Ceo is still not clear. In RefB a small peak at 28 eV was observed and identified 
as the / = 1 cr-plasmon. 

The main goal of this work is to calculate the cross-section of fast electrons scattered by a Ceo molecule and to 
show that the contribution of the plasmons with I = 2 to the cross-section is enough in order to be observed in EELS 
experiments. In the model we use it is possible to separate contributions of different plasmons, i.e. one can obtain 
relative cross-sections. Using this approach we can not calculate plasmon damping. Widths of plasmon peaks will be 
taken from experimental data. 



II. QUANTIZATION OF PLASMA OSCILLATIONS ON A SPHERE IN THE TWO-FLUID MODEL 

We shall consider electrons in the Cgo molecule to be distributed over the surface of a uniformly charged sphere 
— the positive molecular core in the model. The Fourier transformation of the charge density 5n{u>) induced by the 
potential $ including both the external and induced potentials is given by the expression 

a 

where n is the surface tt- and cr-electron density, f a and w„ are respectively the oscillator strength and electron energy 
levels. The index a numerates the electron states. A|| denotes the angular part of the Laplas operator: in the model 
we use electrons cannot move perpendicularly to the surface. 

As the 7r-electron subsystem of graphite is semimetal, we can use the following approximation for the 7r-electron 
contribution to the induced charge density: 



We suggested uo a in equation (^) to be equal to zero and used the sum rule J2 a fair — 1/4 for the oscillator strength 
of 7r-electrons in Cgo- N = 60 and R = 0.4 nm are the number of 7r-electrons and the effective radius of the sphere 
respectively. _ 

On the other hand, the cr-electron subsystem of graphite is dielectric with the gap of order 12 eV (see Rcf.EJ). 
The oscillator strength of the cr-electron subsystem has a maximum at oj a ~ 16 eV. Hence we can use the following 
approximate expression for the cr-electron contribution 



3e 2 N 1 , 
AnmR 2 uj 2 - uj a 2 



Here we used the sum rule J2 a f°t° ~ 3 /4 for the cr-electron subsystem. uj a = 16 eV is the characteristic transition 
frequency of cr-electrons. Is 2 - electrons of carbon atoms are not taken into account in the model. 

The $-field is connected with the total induced charge density 5a = Sn„ + bn a through the Gauss theorem 

-fr Uh+o - -fr Uh-o = 47rS<T - ( 4 ) 
It is convenient to expand the density fields in Fourier series on the sphere 

fry = ai m Yi m ((f>, 6), Sri,, =^2bi m Yi m ((f>,6), a = ^si m Y lm (4),0). (5) 

Ira Ira Ira 

Then using the multipole expansion of the potential (f> 



$ = 



Y,im<t>i™(r/R) l Y lm {cl>,e) : r<R, 
E/ m MW +1 W<M) = r>R 



we can rewrite in the following form 

AttR AttR 



21 + 1 21 + 1 

which gives for (0) and (0) (see also RefsSI^) 



Sim = , 1 ( fl ''» + ( 6 ) 
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a lm = ~fli 2 (ai m + 3bi m ) (7) 

h m = —3Ql 2 (ai m + bi rn ) — UJ a 2 bi m (8) 



where 



a 47reVo Z(Z + 1) TV 
Qi = mi? 2TTT' CT0= 4^- (9) 

It can be easily verified by straightforward calculation that equations (Q) and (JsJ) are the equations of motion 
corresponding to the Hamiltonian 

H = ^ ° Aim,Ai,-m H ^~ - ^ Bi m Bi_ m + 

Ira Im 
Im Im 



111(7 q 



2 

Im 



m Bl- m + Al- m Bl m ) (10) 



where the amplitudes j4; m and -B/ m are the (Z, m)-components of the displacement fields ri and r 2 of the electron 
liquids which arc connected with ai m and 6; m by the following equations 



ao ^ljl+T) 3a y/Z(Z + l) D ri1 , 
a Zrn — "J ^2 ^im, fym — ~ mm- (J-1J 

The coefficients in equation ( |l0| ) are chosen so as to keep it equal to the hydrodynamic Hamiltonian 



^ _ mcr 



J rfdS + ^ J rldS + 3 -^fl J rjdS + I J tSadS (12) 



Hamiltonian ( |Io| ) is nondiagonal in the basis consisting of the functions Yi m (cj), 9) because of the last term. On 
the other hand, this term does not mix amplitudes with different I and \m\ because of the spherical simmetry of the 
model, hence the above Hamiltonian can be diagonalized by the diagonalization of each "block" with a number Im. 

We have shown that the problem can be reduced to the diagonalization of the Hamiltonian 

TOCTq ; 2 , m<J A2 , m<J n 2 A 2 , 3mfJ for, 2 ,2X02,0 O 2 A o /-, o\ 

—^-A H —B H — Sij A H — (3S2; + w CT )B + 2>m<j lli AB (13) 

which can be done by the following transformation 

y/ma A — Qi cos 4>i - Q 2 sin 0; (14) 



y/SmaoB — Qi sin 0/ + Q 2 cos 0/ (15) 

where the angle cf>i can be found from the condition of the disappearance of the nondiagonal term in Hamiltonian jl^ ) 
expressed in terms of the new coordinates Qi and Q 2 - The substitution of ([TJ]) and (|TJ|) into equation ( |l3| ) gives 

2»; 2 +^ 2 

cot2^=-^=^. (16) 
The eigenfrequencies Qij are given by the following equation 

(tti 2 k2 = ~ + a£ ± + W 4 + l 6fi f ) . (17) 

One can clearly see that charge oscillations can be classified by the numbers I and m and the index j corresponding 
to the lower and upper roots of equation (|l7]). The fact that the eigenfrequencies do not depend on the number m is 
due to the spherical symmetry of the considered model. 
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Using (^lj) we can find quantum mechanical operators for amplitudes A\ m and Bi rn . Substituting them into the 
equation for classical charge density fluctuations we obtain the operator corresponding to this quantity 

lm 

where $im,j IS the annihilation operator of the plasmon with numbers lm,j. and the functions g%(4>i) and g 2 (<pi) are 
given by the equations 

gi {<j>i) = cos 4>i - sin <j> u g 2 {<l>i) = cos 4>i + V3 sin (j> h (19) 
the angle cj>i is chosen to be the minimal positive root of equation ( |ici| ) . 



III. FAST ELECTRON SCATTERING 

In this section we will obtain the expression for the scattering amplitude of fast electrons due to charge density 
fluctuations in Ceo- lu the secondary quantization representation the electron-plasmon interaction can be written as 
follows 

H int = [&(-e&)$dr. (20) 

Here * = ^ k e^a^/W is the electron field operator, <l is the potential corresponding to the electron surface density 

distribution Sa((f>, 9), i.e. $ is the solution of the Poisson equation with the charge density S(r — R)8cr(<fr, 9). Hence $ 
and 5& are connected by the following equation 

J |r-r'| 1 ' 

Substituting equation (^l|) and the explicit expressions for the electron field and the electron density fluctuation 
operators (see (|l8|) ) into equaiton (p0|) we obtain the expression for the electron-plasmon interaction 



sr t ( 4 ^) 2e2 /foo 

H in t = - a k'Ok- 7^2 



,k-k'|VV2m 



.iJl(f+i) gjim{Xk _ k > lR)Ylm (|^) (Sim + Sl_J. ! 22 ) 



, V fyj 



We shall consider only the case of fast scattering electrons, i.e. the only relevant processes are those in connection 
with the creation of one plasmon. Hence we can use the Born approximation. The quantitative criterium of the 
applicability of the Born approximation will be discussed further. 

In the initial state all the plasmon occupation numbers are equal to zero. Let us suppose that after scattering only 
the occupation number of the state lm, j is equal to one. Defining the initial and final momenta of the scattering 
electron k and k' respectively we obtain the following expression for the modulus squared of the matrix element of 

Hint 

Kk;0|tf irit |k';*m,j>| 2 = ^^ l M^ g ^ m ^ {qR) Y lm (S) Y L (S) (23) 

where q = k' k. 

As each eigenfrequency f2;j is 21 + 1-fold degenerated, we shall calculate the total cross-section for the transferred 
energy hfl^j. Summing up the contributions of plasmons with given I and j over m we get 

^= / (^| <k; o|^|^im > ^|>(av+n l -a ai )^^^. (24) 
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In the above equation = hk 2 /2m, v is the initial velocity of the electron. 
Using the equation 



J2 Ylm {<\/q)YL (q/<?) = (2/ + 1)/4tt (25) 



I Ti l 



( see Ref.E3) we get after some algebra 



(/ ' T '-< _^' v "„, „,„ , ^ilM^T^^ir 1 k'-k|/?i 



Stt^^ + 1)(2Z + l)«L^Wjfc2 _ —"'-J ^ ^ (26) 



dfi k ' " S, 2 -y H/j V ft Ik'-k' 4 



IV. RESULTS AND DISCUSSION 

In experiments one can measure the number of electrons reaching a detector, i.e. one measures the function J26| ) 
integrated over a certain solid angle. Let this solid angle be restricted by two cones of angles 8\ and 8 2 respectively. 
In this case we can write for the total cross-section <Ji,j{6\,6 2 ) 

|k-l4|fl 

^iA^6 2 ) = m -^l(l+l)(2l+l)g](4> l ) J J -^dx (27) 

Ik-k'Jfl 

where E = mv 2 /2, |k-k<| = \/ k 2 + k' 2 — 2kk' cos8 l , i = 1,2. 

We shall consider the case of small transferred energy, i.e. — <Wk. In principle, several different situations 
are possible. If both the angles 0\ and 8 2 are small, the integration limits in the integral in (p7|) are also small and 
we can expand the spherical Bessel function into a Taylor series and take only the first term of the expansion. As 
j t (x) x l /(2l + 1)!! if x -> 0, we get 

\k-k' 2 \R 

<t u (9i,0 3 )~ J x 2l ~ 3 dx (28) 

Ik-k'JB 

and one can clearly see that the contribution of the plasmons with I = 1 is dominant. In the opposite situation 
|k — k^|i? 3> 1, i = 1,2 we can use the asymptotic expression for the spherical Bessel function ji(x) = sin (x — irl/2)/x 
if x — > oo. As the integration limits do not depend on I significantly, we can see that the values of the integral in ( 571 ) 
are of the same order for different I. Hence the relative contributions of plasmons with different I depend only on the 
factor gf A<f>i) In the case |k — k^|iJ ~ 1, i = 1,2 it is possible to estimate the integral in ( |27| ) by the substitution 
ji(x) ~ x I (21 + 1)!!. One can clearly see that only the contribution of the plasmons with small I is relevant because 

of the rapidly increasing multiplier ((21 + 1)!!) in the denominator. 

Strictly speaking the above analysis is not applicable for the plasmons with I > 3. The Ceo symmetry group is 
icosahedral and it does not have irred ucib le representations with a dimension of more than five. Hence the plasmons 
with I > 3 in the spherical model splitE2rt3 and the classification of plasmons in the real Ceo molecule is different. We 
will not discuss further this difference because in the experiments which we will analyse only the contribution of the 
plasmons I = 1, 2 is relevant. On the othep-hand, the small size of Ceo leads to the strong damping of plasmons. The 
estimation for the spherical model in Ref.E3 shows that the ulasmons with I > 3 could hardly exist. 

We shall compare our results to the experiment from Ref.El . In this work the number of electrons with E = 1 keV 
scattering in the solid angle 0.5° < 8 < 2.5° (i.e. 9% = 0.5°, 8 2 = 2.5°) per unit time was measured. The average 
cross-section 

= <J Lj (8 1 ,8 2 )/27r(cos8 1 ~ cos# 2 ) (29) 

obtained using the spherical model for the plasmons with I = 1 ... 5 is shown in Table |[ We used the value u> a = 16 eV 
for the characteristic cr-electron transition frequencies and R = 0.4 nm for the radius of the sphere in the model. For 
the conditions of the analised experiments |k — k^|i? ~ 1, i — 1,2 and one can clearly see that the contribution of the 
plasmons with I > 2 decays rapidly. This is in good agreement with the qualitative consideration. The I = 1 ir- and 
CT-plasmons and the quadrupole cr-plasmon make the most significant contribution to the cross-section. 
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Plasmons in Cqo are characterized by strong damping due to the small size of the molecule: in small objects 
the energy of plasma oscillations per one electron hft/N is strongly enhanced and leads to strong electron-electron 
scattering, meanwhile for bulk plasmons this energy is negligible. As energies of ir- and cr-plasmons do not depend on 
I significantly (see Table|j), the parameter UCl/N is approximately the same for n- and cr-plasmons respectively and 
we can expect that the damping of plasmons does not depend on I strongly. 

In the above calculation of the cross-section we did not take into account plasmon damping. The characteristic 
inverse life time rjj taken from the experiment is equal to 1.5 and 8 eV for 7r- and cr-plasmons respectively. The 
function 



is plotted in Fig. 1 together with the experimental curve. One can see that our results are in good agreement with 
the experiment. 

As it can be seen from Table Q and Fig. 1 the first peak on the graph corresponds to the dipole 7r-plasmon; the 
second peak corresponds to the dipole and quadrupole cr-plasmons. The contributions of the other plasmons are small. 
It has already been discussed that ratios 0^1/0^2 are proportional in our case to ratios (gf 1 /n^i)/(gf 2 /rti,2)- Note 
that although fij 1 < fi;^ we have for the cross-sections 07,1 > 07,2 and the ratio o^i/cr/^ decreases if I increases. This 
is the consequence of the fact that gi.i — > if I — ► 00 : for large I we can write using equation jl6| ) cot (f>i — ► 
hence 137,1 — > (see (|Ts| ) ) and we obtain an interesting result: the ratio 07.1/07.2 tends to zero if I — > 00. 

Let us discuss the.-c.ondition of the applicability of the Born approximation. The Born approximation is valid if 
\U\a <C hv (see Ref.E3) where \U\ and a are the characteristic value of the potential and its radius respectively, v 
is the velocity of the scattering particle. In our case the characteristic value of charge density fluctuations 5a is 
(e/ a 2 ) y/Ti(T 1 mil where O ^ 10 eV is the characteristic oscillation frequency. Hence \U\ ~ (e 2 / a) WUoq / mil and the 
condition of the applicability can be written as follows 



Substituting ctq = N/AuR 2 , a ~ i?, R = 0.4 nm into ( [31] ) we see that the Born approximation is applicable if the 
velocity of the scattering electron satisfies the condition v 3> 2 • 10 8 cm/s. In the analysed experiments the electron 
velocity is 1.88 • 10 9 cm/s, i.e. the Born approximation is applicable. !_. 

Briefly summarizing we have shown that the simple two-fluid model used by the authors oi!3 for the predictions of 
the collective spectrum of Ceo turned out to be good for the quantitative calculations of the intensities of the plasmon 
peaks in EELS spectra. The advantage of the description based on this model is that it is physically and technically 
simple and the model itself does not include any fitting parameters (with the exception of the decay rates which 
however can not be properly calculated in more sophisticated techniqueaj) . All this makes the model useful for the 
quantativc consideration of the collective behaviour of electrons in more complicated fullerene systems. 

TABLE I. Frequencies and intencities of plasmon peaks calculated using formulae (17) and (27) for IkeV incindent electrons 




(30) 




(31) 



at 1.5° ± 1°. 



Type of plasmon 



Plasmon frequency, eV 



Cross section, nm 2 /sr 



1=1, j=l 
1=1, j=2 
1=2, j=l 
1=2, j=2 
1=3, j=l 
1=3, j=2 
1=4, j=l 
1=4, j=2 
1=5, j=l 
1=5, j=2 



5.95 
22.21 
6.68 
26.49 
7.03 
30.12 
7.23 
33.35 
7.36 
36.29 



8.4 • 10~ 3 

0.22 
5.6 ■ 10~ 4 

0.024 



3.82 
9.56 
0.69 
5.16 
0.09 
1.35 



G 



1.2 
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